A square lattice algebraic spin liquid with SO(5) symmetry 
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We propose a critical spin liquid ground state for S = 1/2 antiferromagnets on the square lat- 
tice. In a renormalization group analysis of the 'staggered flux' algebraic spin liquid, we examine 
perturbations, present in the antiferromagnet, which break its global SU(4) symmetry to SO(5). At 
physical parameter values, we find an instability towards a fixed point with SO(5) symmetry. We 
discuss the possibility that this fixed point describes a transition between the Neel and valence bond 
solid states, and the relationship to the S0(5) non-linear sigma model of Tanaka and Hu. 



Critical spin liquids appear in a variety of contexts 
in studies of correlated electrons in insulators and su- 
perconductors. These are states in vi^hich spin rotation 
symmetry is preserved and there is a gapless spectrum 
of spin excitations which do not have a quasiparticle in- 
terpretation. The simplest, and best understood, states 
[l|, Q are described by the Wilson-Fisher fixed point of 
the Landau- Ginzburg- Wilson theory of fluctuations of 
the vector antiferromagnetic (Neel) order parameter N . 
These provide a description of the quantum phase transi- 
tion of dimerized antiferromagnets between a state with 
long-range Neel order and a spin-gapped state with S = 1 
quasiparticle excitations which are quanta of the field TV. 
However, the anomalous dimension of the field N is quite 
small at the critical fixed point, implying that a pertur- 
bative description in terms of TV quasiparticles provides a 
reasonable description of the zero temperarure spectrum. 

A separate category of critical (or 'algebraic') spin liq- 
uids involve a description in terms of neutral S — 1/2 
excitations ('spinous'). The first example of this was the 
'staggered flux' spin liquid state of Affleck and Marston 
whose low energy theory involved spinous, repre- 
sented two-component massless Dirac fermions (a = 
1 ... TV = 4 is a combined SU(2) spin and valley index), 
coupled to an emergent U(l) gauge field a^. The spin 
fluctuations of this theory have been better understood 
in subsequent work 0, [^, 0]: as we will review below, for 
TV sufficiently large, the low energy theory is a confor- 
mal field theory (CFT) with a global SU(4) symmetry. 
Another spin liquid state involves bosonic spinous rep- 
resented by relativistic scalars (p = 1,2 is a SU(2) 
spin index) coupled to a U(l) gauge field [8]. It has been 
proposed [9] that this is realized as a CFT describing 
a quantum critical point between states with Neel and 
valence bond solid (VBS) order. 

This paper will propose a critical spin liquid ground 
state described by a CFT with a global S0(5) symmetry. 
We begin with the SU(4) CFT of Dirac fermions noted 
above, and examine the renormalization group (RG) fiow 
of all perturbations which preserve relativistic invariance 
and at least a global S0(5) C SU(4) symmetry. The 
global symmetry of the underlying antiferromagnet in- 
volves only a continuous SU(2) spin rotation symmetry 
and various discrete space group symmetries, and all such 



S0(5) invariant perturbations will generically be present. 
Using an expansion defined below, we find that for phys- 
ical parameter values the SU(4) fixed point is unstable 
to fiow towards a S0(5) invariant fixed point. 

Our motivation for examining CFTs with SO (5) sym- 
metry comes from an interesting proposal by Tanaka 
and Hu [l^ (see also the work of Senthil and Fisher 
[lH). They examined the quantum fluctuations of the 
3-component Neel order parameter, TV, and the com- 
plex VBS order parameter S, and suggested that they 
be combined into a single five-component real vector E^, 
with a = 1 ... 5 which transforms as the fundamental of 
an enlarged SO (5) group. In a spin liquid with such an 
SO (5) symmetry, the anomalous dimensions of TV and 2 
would be equal, and Sandvik's numerical results [l3| on 
the quantum critical point between the Neel and VBS 
states are consistent with such an equality. 

We begin by reviewing the SU(4)-invariant CFT of 
Dirac fermions, largely following the notation of Ref. 0- 
The CFT is described by the Euclidean spacetime action 
So — J (PrdrCo, where 

/:o=^V(5p+m^)«'a (1) 

where fi = t, x,y is a spacetime index, 7^ are the Dirac 
matrices, 4* = 4*^70 and is an emergent U(l) gauge 
field. As shown in earlier work 0, 0], this action de- 
fines an SU(TV) invariant CFT in an expansion in 1/TV. 
The combined Neel- VBS operator, E^, can be written in 
terms of the by 

Ea = ^Fa* (2) 

where F^ are five 4x4 matrices from the SU(4) al- 
gebra. This algebra can be realized using the ten- 
sor product of two independent sets of Pauli matri- 
ces, /X and <T, and Hermele et al. showed that Fq = 
{fi^a'^ , li^a^ , fi^a^ , fi'^ , IJ,^). A curious, and key, property 
of the Fa is that they anti-commute, {Fa, Ff,} = 2Sab, and 
so they are Dirac matrices of five spacetime dimensions. 
The 10 generators of the SO (5) group, under which Ea 
transforms as a SO (5) fundamental, are obtained from 
the commutators of the Fa: 

Fa6 = ;^[ra,Fb]. (3) 



2 



The Ta and Tab are the complete set of SU(4) generators. 

It will be important for our analysis to be able to gener- 
alize these order parameters, and the associated algebraic 
structure, from SU(4) to general SU{A^), so as to allow a 
systematic expansion. A similar strategy was used 
in the context of chiral symmetry breaking of three di- 
mensional QED [is'l. However, the above embedding of 
S0(5) into SU(4) relies on the spinor representations of 
SO (5), and this does not have a suitable generalization. 
However, we note that there is an antisymmetric matrix 
J = iay^'^, with J'^ = —1, under which 



J^abJ — 



(4) 



for all ah. Eq. (|4|) is the defining relation for generators 
of the Sp(4) subgroup of SU(4), and we have just estab- 
lished the well-known congruence Sp(4) = SO (5). The 
embedding of Sp(A) into SU(A) generalizes easily to all 
even N , with an N x N antisymmetric matrix obeying 
J"^ = -I. We will therefore study here the SU(A) in- 
variant CFT in Eq. ([T]) with a = 1 ... A, while allowing 
perturbations which are invariant under Sp(A). 

A linear stability analysis of this SU(A) CFT has been 
carried out earlier |7| for a hmited set of perturbations. 
For sufficiently large N, all perturbations are believed to 
be irrelevant. However, the anomalous dimensions aris- 
ing at order 1/A can be quite large, and we shall show 
find below a perturbation which becomes relevant when 
its scaling dimension is evaluated at A = 4. We are 
also interested in finding a systematic approach to deter- 
mining the fate of such a relevant perturbation, beyond 
a linear stability analysis. To this end, we will allow 
the tree-level scaling dimensions to vary as a function of 
spatial dimensionality, d, as is common in other critical 
phenomena contexts. With Dirac fermions, there is the 
subtle issue of dimensional continuation of the Dirac ma- 
trices, 7^*; as is commonly done we will deal with 
this by applying the Dirac algebra and phase space fac- 
tors as in d = 2. Our stablility analysis of spin liquids 
and their perturbations is formally justified by taking 
(d — 1) oc 1/A, and then expanding in 1/A. 

It is also interesting to consider application of this 
method to antiferromagnets in d = 1. Although we 
will not describe the computation here, it is necessary 
to adapt our results to Dirac matrices in c? = 1. From 
such a computation, we reproduced the results of Affleck 
[l3 | on the spectrum of scaling dimensions of operators 
with SU(A) and Sp(A) symmetry at the fixed points 
described by WNZW models. 

We now present our RG results for perturbations of the 
CFT in Eq. ([T]). We begin by considering perturbations 
which are invariant under SU(A). To the order we are 
working, there are only two independent perturbations, 
which we write as 



^(vl/"7'^vl/„)(vl/^^*^) (5) 



FIG. 1: Feyriman diagrams contribute to the linear orders 
in both Eqs. |6] and [O] The dashed lines are dressed photon 
propagators, and the full circles denote the trace in Dirac 
space. 





FIG. 2: Feynman diagrams which only contribute to the linear 
orders in Eqs. [H but not in Eqs. El 



where the circular brackets indicate a trace over in- 
dices in Dirac space. Other possible terms, such as 

a a 

(vp 'i'0){^^ \E'q) and 7^\1/^)(\E' 7^*q), can be shown 
to be linearly related to the terms in Eq. ([5|). 

From the diagrams shown in Fig. [U Fig. [5] and Fig. [31 
we obtained the following RG equations for a rescaling 
by a factor e^: 
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In the terms linear in the A on the right-hand-side, we 
have computed co-efficients to order 1/A, and the 1/A 
corrections come from the dressed photon propagator Q 
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FIG. 3: Feynman diagrams which contribute to the quadratic 
order of the RG equations (O and Notice that since we 
only calculate to the order of unity in the quadratic terms, 
diagram G only contributes to equation ((6} but not (|9]), and 
diagram H only contributes to equation @ but not ((SJ. 



For the terms quadratic in A to be of the same order as 
the linear terms, we need only compute the co-efScients 
to order unity, as is the case above. 

The RG equations in Eq. ^ have several fixed points, 
but we begin by considering the fixed point at Ai = 
A2 — 0. The eigenvalues at this fixed point are 1 — d — 
(128 ± 64V7)/(3A^7r2). At the physical values of d = 2 
and N — 4, these eigenvalues evaluate to —0.651 and 
—3.510. So both are negative and the Ai = A2 = 
fixed point is stable. None of the other fixed points of 
Eq. ^ were found to be stable at these values of d and 
N. By examining the A'^ dependence of the eigenvalues 
at Ai = A2 = we conclude that the SU(7V) CFT defined 
by Eq. ([T]) is stable to SU(-/V)-invariant perturbations for 
N > 1.40/(d- 1). 

Next we consider the additional perturbations of £0 
when the global symmetry is reduced from SV{N) to 
Sp(A'^). A simple analysis shows there is only one allowed 
term 



C2 = ^Ja^J^\^"^f3)i^''^5). 



(8) 



A second possible term Ja'yJ^^i'^"j'^^p)i'^^jf^^5) re- 
duces to the above term after application of Ficrz iden- 
tities. 

From the diagrams in Fig. [1] and diagram H in Fig. [3l 
the RG equations for £2 reads (Notice that diagrams in 
Fig. [2] and diagram G in Fig. [3] do not contribute to the 
leading order of 1/N expansion) 
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This has fixed points at A3 = and A3 = A3 = 37r (1 — 
d + 6i/{Nn^)). At d = 2 and iV = 4 we now find a 
result which is very different from the SU(A^) perturba- 
tions above. The A3 = fixed point is unstable with RG 
eigenvalue 0.621, while the fixed point at A3 = A| > 
is stable; for general N, we find that the stablity of the 
A3 = A^ fixed point holds for N < 6.48/(d- 1). So for 
1.40/(d-l) <N < 6.48/(d-l), the theory Ca+Ci+C2 
flows to a fixed point with Ai = A2 = and A3 = AJ 
which describes our advertised Sp(iV)-invariant critical 
spin liquid. 



FIG. 4: Feynman diagrams which contribute to the difference 
of scaling dimensions of fermion bilinears ^'Fa*!' and v&Fai,*!'. 



The scaling dimensions of all 16 fermion bilinears 
fTa* (Ta are SU(4) generators with a = 1, • • • 15) equal 
at the fixed point with Ai = which respects the SU(4) 
symmetry. At the order of 1/N, the scaling dimensions 
read: 



A(^'ra*) = 2 - 
A(f*) = 2 
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with = 4, the difference between the two scaling di- 
mensions above is from the diagrams similar to the ones 
in Fig. [3] with two photon propagators and a trace 
in the fermion flavor space, which only contributes to 
fermion bilinear 5'^'. At the Sp(4) symmetric fixed point, 
the scaling dimensions of fermion bilinears are classified 
as the_representation of Sp(4)~SO(5) group: ^Ta^ 
and ^'Fafc^' form scalar, vector and adjoint representa- 
tions of S0(5) group respectively, and the scaling dimen- 
sions of fermion bilinears within the same representation 
are equal to each other. 

For larger A'^, the scaling dimensions of the fermion bi- 
linears at the Sp(A'^) fixed point deviate from their value 
at the SU(7V) fixed point at the order of and re- 

quires a lot more calculations. But their differences at 
1/A^^ order can be calculated readily from diagrams in 
Fig. HI 



A(4'r,*) - A(«'r,fc*) 
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(11) 



Here F^ and Tab together form a fundamental represen- 
tation of SU(iV) algebra, and Tab form the spinor repre- 
sentation of Sp(A^) subalgebra. 

To fully analyze the physical implications of this fixed 
point, we have to examine the fate of all perturbations 
which further reduce the global symmetry from Sp(4) 
down to those required by the SU(2) spin rotation sym- 
metry and the square lattice space group. There are a 
large number of such additional perturbations, and an- 
alyzing them all would require an analysis of daunting 
complexity. We also need a procedure for generalizing 
such perturbations to general Sp(iV) operators to enable 
a 1/N expansion, and there is no unique and natural 
choice like the one we have used so far; the results will 
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depend upon the particular choices made for the invari- 
ant subgroups of Sp(A^). We will therefore not present 
such an analysis here. Additional perturbations which 
break Lorentz invariance are also possible; there were ex- 
amined by Hermele et al. and found to be irrelevant. 

Should no relevant perturbations emerge at the Sp(4) 
fixed point, it would describe a stable critical spin liquid 
phase. Otherwise it would be a (multi-) critical point 
between ordered phases, with the dimensionality of the 
phase diagram determined by the number of relevant op- 
erators. An intriguing possibility is that there is only one 
relevant perturbation, which drives the system to a Neel 
or a VBS state on opposite sides of the Sp(4)-invariant 
critical point. 

Such a Sp(4) = SO (5) fixed point separatin g N eel and 
VBS states was suggested by Tanaka and Hu [10|. They 
further proposed a SO (5) non- linear sigma model, with 
a Wess-Zumino term which could realize that a critical 
state. However, our Sp(4) critical point also has a U(l) 
gauge field, and an associated conserved topological cur- 
rent, and there is no analog of this conserved current in 
the Tanaka- Hu sigma model. So it is likely that our Sp(4) 
critical spin liquid is distinct from their proposal ^15j . 

A large number of possible spin liquid ground states 
have been proposed for the square lattice antiferromag- 
net. All previous proposals have been associated with a 
mean-field saddle point of a theory of electrically neutral 
spinous which are either fermions or bosons. This paper 
has proposed a novel type of a spin liquid, which does not 
have a direct mean-field realization, but is induced by the 
gauge fluctuations about a mean-field saddle point. The 
only numerical evidence so far of a spin liquid state on 
the square lattice for SU(2) antiferromagnets is in the 
studies of the transition point between Neel and VBS 
Our SO (5) spin liquid is a candidate 



observable operators 0| defined by "^Tab^ all have equal 
scaling dimensions, which are distinct from those of the 
Neel and VBS orders. 

We thank T. Senthil for valuable discussions. This 
work was supported by NSF Grant No. DMR-0537077. 
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for this state, as it can explain the possible equality of 
the scaling dimensions of the Neel and VBS operators. A 
further testable property of our spin liquid is that the 10 
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